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1. Introduction 

Attempts to tackle to isotropization problem in relativistic cosmology date at least since 
the pioneering work of Collins and Hawking [Q who proved isotropization theorems for 
certain classes of orthogonal Bianchi spacetimes. A further advance is the recent work 
by Heusler [0] which extended the Collins-Hawking results to the case of a minimally 
coupled scalar field matter source with a general convex potential. Interest in a 
particular approach to the isotropization problem in cosmology renewed after the advent 
of inflation as a mechanism for solving the problem and focused on proving the so- 
called cosmic no-hair conjecture. This conjecture roughly speaking states that general 
cosmological initial data sets when evolved through the gravitational field equations are 
attracted (in a sense that can be made precise) by the de Sitter space of inflation. In 
other words, if this conjecture is true, inflation is a 'transient' attractor of such sets, 
these in turn quickly isotropize after the inflationary period and the inflationary regime 
could thus be regarded as 'natural' in view of its prediction of the (observed) large-scale 
homogeneity and isotropy of the universe. 

Higher derivative quantum corrections to the gravitational action of classical general 
relativity are generally expected to play a significant role at very high energies where 
a quantum gravitational field will presumably dominate. It is not unreasonable to 
consider classical cosmology in theories coming out of such nonlinear (higher derivative) 
gravitational Lagrangians and in fact, one expects that there exist close links between 
properties of such "higher derivative cosmologies" and those of general relativistic 
cosmology. It is obvious that the resolution of aspects of the singularity, isotropization 
and recollapse problems of cosmological spacetimes is of paramount importance also in 
this extended framework. 

In this paper, we prove the cosmic no-hair conjecture for all Bianchi (including type 
IX) matter cosmologies which are derived from a quadratic action in the scalar curvature 
that is of the form L = R + (5B? + Lm where Lm denotes the matter Lagrangian. 

Other aspects of this problem were considered previously, for instance, for the case 
of the Einstein field equations with a cosmological constant by Wald , and Maeda [^ 
for the case of the quadratic theory L = R+(3B? in vacuum for all Bianchi models except 
the IX. More specifically, Wald proved, in the particular case of a true cosmological 
constant, that this conjecture is true for all Bianchi models, except possibly for type IX 
. Wald's proof served as the prototype for most subsequent works where, either a true 
cosmological constant was assumed or the validity of the cosmic no-hair conjecture 
was tested in the context of particular models: chaotic inflation |^, , power-law inflation 
0, H, inflation due to a E? theory [^, |^, |1^ or in more general / (i?) theories [|lT . 



Our results contribute to an extension of previous works in two ways. We provide 
a detailed study of the role that matter plays in the problem and treat rigorously the 



important case of the Bianchi IX model (including matter fields) in the case of the full 
self -interacting potential arising from the conformal transformation. Our proof relies 
on the use of an energy-like function which we call the Moss-Sahni function and is 
accomplished through a detailed study of energy estimates of the associated quantities 
that obey certain differential inequalities. We show that once the required bound for the 
Moss-Sahni function is obtained, the spatial part of the scalar curvature asymptotically 
approaches zero almost exponentially. 
The field equations take the form 



nab ^9abn ^^2/3i? 



2VaVbi? - 2gabaR + -QabR' 



Under a conformal transformation of the metric 

^„6 = (1 + 2PR) gab, 
with 



Tab (g) . (1) 



(2) 
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■ In (1 + 2(3R) 



the field equations ([^) become the Einstein equations in the new space-time (M, g) 



R. 



■ab 



-gab R = Va'-P'Vb^P - -gab (VcV^VV) " 9abV + Tab (g) 



(3) 



(4) 



U^ - V' {^) = 0, 
where the potential V is [ p!2| , 



V 



1(3 



1 — exp 




(5) 



(6) 



As Maeda 0] has pointed out, this potential has a long and flat plateau. When ip is 
far from the minimum of the potential, V is almost constant Voo = lim^^+oo V {(p) = 
1/ [8(5) . Thus V has the general properties for inflation to commence and \4o behaves 
as a cosmological term. The constant /? is of order 10^^ /p^ |I3| . 

The plan of the paper is as follows. In the next Section we prove the cosmic no-hair 
conjecture for all Bianchi universes except IX, while Section 3 is devoted to the analysis 
of the Bianchi IX model. Precise statements of the results proved may be found in those 
Sections. We conclude by pointing out how our results could lead to the treatment of 
several more general cases. 



2. All Bianchi models except type IX 

In this Section we show using Einstein's equations (|^) and the equation of motion of 
the scalar field (^), that, all Bianchi models which are initially on the flat plateau of the 



potential (j^, except probably Bianchi IX, with a matter content satisfying the strong and 
dominant energy conditions, rapidly approach de Sitter space-time. This is established 
for the full scalar field potential arising from the conformal transformation. 

We shall work exclusively in the conformal picture and so for simplicity we drop 
the tilde. Choose n°' to be the unit geodesic vector field, normal to the homogeneous 
hypersurfaces. The spatial metric is related to the space-time metric as hah = Qab + ^a^6 
and rf-V a = d/dt, where t denotes proper time along the integral curves of n°. The 
scalar curvature R becomes a function of time, hence the scalar field introduced by eq. 
(§) is homogeneous. Ordinary matter is assumed to satisfy the strong and dominant 
energy conditions, namely, 

Ts = {Tab-lTgab)n''n'>0 

Td = Tabu"'u'^ > 0, Tl^u^ non — spacelike, 

for any unit timelike vector field m°. 

We make use of the time-time component of Einstein's equations (§) 

lK' = a' + Td + l^' + V-l^'^R (8) 

and the Raychaudhuri equation 

K= -\k^ - 2a^ - Ts- ^^ +V. (9) 

The equation of motion for the scalar field (^ becomes 

(p + K^ + V'{^)=Q. (10) 

Here K is the trace of the extrinsic curvature Kab of the homogeneous hypersurfaces 
and is related to the determinant h of the spatial metric by 

ir=|(ln/.V^). (11) 

We note that in P, |l4|, the system above is reduced to a two dimensional autonomous 



form as a consequence of the exponential form of the potential which we do not assume. 
Our case is more involved since the system cannot be reduced further and so we have 
to proceed in a different way. We make the following assumptions: 

1. The initial value ipi of the scalar field is large and positive (i.e., the universe is on 
the fiat plateau of the potential). 

2. The kinetic energy of the field is negligible compared to the potential energy. 

3. The universe is initially expanding i.e., K > at some arbitrary time tj. 

With assumption (3), Eq. (|) implies that the universe will expand for all subsequent 
times i.e., K > {] for t >ti for all Bianchi models, except possibly type- IX. Defining the 
energy density of the scalar field hj E = ^ ^ +V, we find by (|T0|) that 

E = -Kif'^. (12) 



Hence, in an expanding universe the field looses energy and slowly rolls down the 
potential. The "effective" regime of inflation is the phase during the time interval 
tf — ti needed for the scalar field to evolve from its initial value ipi to a smaller value 
ipf, where iff is determined by the condition that V (iff) ^ rjVoo- The numerical factor 
1] is of order say 0.9, but its precise value is irrelevant. As we shall show, the universe 
becomes de Sitter space during the effective regime followed by the usual FRW model 
when the cosmological term vanishes. To this end, we define a function S, which plays 
the same role as in Moss and Sahni [^], by 

S = -K'^ - E. (13) 

In all Bianchi models except IX, using (8) we see that this function is non-negative due 
to the dominant energy condition and the fact that in these models the scalar spatial 
curvature is non-positive. Taking the time derivative of S and using eqs. (||) and (j^) 
we obtain 

S = ~KS-'^{2a^ + Ts)K. (14) 

It follows that S < -{2/3)KS or, 

S < -^S^3{S + E). (15) 

This differential inequality cannot be integrated immediately because E is a function of 
time (albeit slowly- varying) . However as is well known ||16|, in order to have inflation, 
E must be bounded below and in this case assumption (1) above implies that we may 
assume that the scalar fleld is large enough and so, without loss of generality, we may 
set E > rjVoo- The choice of t] affects only slightly the isotropization time, but the 
qualitative behavior of the model is the same. Therefore, inequality ([l5|) implies that 

S < — TT^-T, ^ = \^VVoo/S. (16) 

smh [mt) 

From eq. (^ we see that, as t increases, the shear, three-curvature and energy density 
of matter rapidly approach zero and because of the dominant energy condition, all 
components of the energy-momentum tensor approach zero. It follows that the universe 
isotropizes within one Hubble time (1/a/\4o ~ 10^ ^pl)- 

3. Bianchi type IX models 

In this Section we show that Bianchi-type IX also isotropizes if initially the scalar three- 
curvature ^^^R is less than the potential V of the scalar field. In the present case, the 
scalar curvature ^^^R may be positive and the argument of the previous Section does not 
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apply directly. In particular, the function 5* is not bounded below from zero.f However, 
we can estimate an upper bound for S by observing that, either S is bounded above 
from zero, or, if S is initially positive, inequality ( jIS]) implies as before that an upper 
bound for S is given by (|16|). Collectively, we have the overall result that 



S < raaxlO, 3?7i^ sinh '^ (mt)> . (17) 

(Had we chosen an exponential potential, this bound would have had the form found in 
0, p. 1418). An estimation of a lower bound for S derives from the fact that the largest 
positive value the spatial curvature can achieve is determined by the determinant of the 
three metric 

(=^)i?„,ax oc /i-^/^ = exp (-2a) . (18) 

We obtain a lower bound for 5* in the following way. We assume that the three- 
scalar curvature ^^^R is less than the potential energy since otherwise it is known that 
premature recoUapse commences before vacuum domination drives the universe to inflate 
ir?! . Then, eq. (H) immediately yields 

^K^ _ i <^V 0. 
3 2 



Observe that for any < A < ^2/3 the above inequality implies that for K > 0, 
U. 1, . 1,. 1 



-K --X^> -K -^ 
3 2-3 y/6 



V 



> 0. (19) 



Suppose that initially i.e., at time tj we have V > ^^-"-Rmax- We claim that this inequality 
holds during the whole period of the effective regime of inflation. Following |Q we define 
a function / by 

/(t)^ln-|— , te[U,tf], (20) 

whose initial value is positive by the above assumption. Differentiating we find 

exp ( — 1 



f=^ + 2c- 



I — exp f 




(21) 



since K = 3 a, hj eq. ([TT|) . During the effective regime of infiation V{!f) > rjVoo, 
hence solving for (p this inequality we immediately verify that the coefficient of (f in the 

f However, as Wald [pi points out, for some non-highly positively curved models premature recollapse 
may be avoided provided a large positive cosmological constant compensates the —^^^^R term in eq. 
(^. In our case, of course, it is the potential V (</?) which acts as a 'cosmological term'. A similar proof 
to that of Wald [^ for the Bianchi IX model, was given by Kitada and Maeda |§] in the case of an 
exponential potential leading to power-law inflation. 



brackets is less than 1/\/Q. Therefore, inequahty (p^Qf ) imphes that / > . We conclude 
that, ii K > Q and / (tj) > initially, then for t G [ti,t/] we obtain / (t) > / (tj) and 
our assertion follows. 

We are now in a position to estimate the required bound for S. First, from eq. (H) 
it is evident that 

and since V/2 > (^^_Rmax/2, eq. (H) again implies that 



K>^lv> ^^r^\/oo. (23) 

Remembering that K = 3a and using the last inequality and eq. (0), we see that 
^^-"-Rmax decays faster than exp —J2'r]Voo/'3 {t — ti) . Therefore ( ^2]) becomes 



■^(^)i?_(t.)exp 




< S. (24) 



Combining this lower bound for S with the upper bound (|T7|)) we conclude that S 
vanishes almost exponentially. From eq. (^, —25 < *^'^^i? < ^'^^-Rmax and therefore '^^•'i? 
damps to zero just as S and '-^■'-Rmax- 

We conclude that when the universe is initially on the plateau, the shear, the scalar 
three-curvature and all components of the stress-energy tensor approach zero almost 
exponentially fast with a time constant of order ~ l/VKx^-t 

We discussed inflation in the equivalent space-time (M, g) , but it is not obvious 
that the above attractor property is maintained in the original space-time (M, g) . This 
is probably an unimportant question since there is much evidence that in most relevant 



cases the rescaled metric g is the real physical metric |[15[ . However, Maeda has pointed 
out m that inflation also naturally occurs in the original picture since during inflation 
the scalar field changes very slowly and the two metrics are related by g = exp (v/|</5J g. 
We now move on to show that the time needed for the potential energy to reach 
its minimum is much larger than the time of isotropization. Therefore, the universe 
reaches the potential minimum {(f = 0) at which the cosmological term vanishes and 
consequently evolves according to the standard Friedmann cosmology. Wald's proof 
assumes the existence of a cosmological constant. However, in realistic inflationary 
models the universe does not have a true cosmological constant but rather a vacuum 
energy density which during the slow evolution of the scalar field behaves like a 
cosmological term which eventually vaniches. Therefore we are faced with the question 
of whether or not the universe evolves towards a de Sitter type state before the potential 

f Actually, the time constant for isotropization in type IX is longer by v2 than in other types. The 
situation is similar to that encountered in Kitada and Maeda M . 



energy of the scalar field reaches its minimum. In any consistent no-hair theorem one 
has to verify that the time necessary for isotropization is small compared to the time the 
field reaches the minimum of the potential. The following argument shows that in our 
case the vacuum energy is not exhausted before the universe is completely isotropized. 
(Similar estimates were derived in for the case of the standard quadratic potential.) 
Imagine that at the beginning of inflation the universe is on the fiat plateau of the 
potential. It is evident that the time tf — ti needed for the scalar field to evolve from its 
initial value (fi to a smaller value (pf is smaller in the absence of damping than that in 
its presence. In the absence of damping, tf — ti is easily obtained from integrating the 
equation of motion of the scalar field ( p^Of ) {K = 0). Taking for example a iff such that 
V {'■Pf) — rjVoo one finds that tf — ti is more than 65 times the time r of isotropization. 
The presence of damping increases the time interval tf — ti and a larger anisotropy 
damps more efficiently the slow rolling of the scalar field, thus producing more infiation. 
It follows that when due account of damping is taken, the period of the effective regime 
of inflation is more than sufficient for the complete isotropization of the universe. 

4. Conclusions 

It is known that the quasi-exponential solution of the R + [3R^ theory considered 
in this paper is an attractor to all isotropic solutions of this theory. Our analysis here 
suggests that the solution could be shown to attract the class of orthogonal Bianchi 
universes as well. If this is indeed the case, it would be interesting to ask whether this 
attractor is unique in the space of all higher order gravity theories. One way to directly 
verify this might be through a perturbation analysis and the corresponding investigation 
of the asymptotic structure of the solutions. 

We believe that the proof of the cosmic no-hair conjecture presented here could 
be extended with adjustments in two directions namely, in an arbitrary number of 
dimensions and also to the class of titled Bianchi models. In fact, such a demonstration 
in some titled models would amount to a first test of this conjecture in cases of some 
inhomogeneity. An analysis along these lines might be more tractable than say attacking 
directly a genuine inhomogeneous case such as for instance that of G2 cosmologies 
wherein the dynamics is described by systems of partial differential equations. Also 
a more direct analysis along the lines adopted here could lead to an extension of our 
results in higher dimensions. 
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